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1 INTRODUCTION 



Exactly solving the MLLA evolution equations for the quark and gluon inclusive spectra and for 2- 
particle correlations inside one jet provided, at small x, in [1], analytical expressions for these ob- 
servables, which were unfortunately limited, for technical reasons to the "limiting spectrum" A = 
ln(Qo/ Aqcd = 0. The goal of this second work is to go beyond this limit in an approximate scheme 
which proves very economical and powerful: the steepest descent (SD) method. It offers sizable technical 
progress in the calculation of both observables. 

First, we perform a SD evaluation of the (quark and) gluon single inclusive distributions. Their full 
dependence on A is given, including the normalization. The well known shift to smaller values of x of 
the maximum of the distribution, as compared with DLA calculations is checked, as well as its Gaussian 
shape around the maximum. Comparison with the results obtained numerically in is done. 

As shown in knowing the logarithmic derivatives of the inclusive spectra immediately gives access 
to 2-particle correlations. This is accordingly our next step. Since, in particular, the former prove to be 
infra-red stable in the limit A — > 0, the result can be safely compared with the exact one obtained in 
The agreement turns out to be excellent, and increases with the energy scale of the process. 

Last, we evaluate 2-particle correlations inside one high energy jet and study their behavior at Qq ^ 
Aqcd- That one recovers the results of Fong & Webber [3| close to the peak of the single inclusive 
distribution and when A — > is an important test of the validity and efficiency of the SD method. The 
quantitative predictions do not substantially differ from the ones of [ 1 1 for the "limiting spectrum", which 
stays the best candidate to reproduce experimental results. 

A conclusion summarizes the achievements, limitations and expectations of [ 1 ] and of the present work. 
It is completed with two technical appendices. 

2 STEEPEST DESCENT EVALUATION OF THE SINGLE INCLUSIVE 
DISTRIBUTION 

We consider the production of one hadron inside a quark or a gluon jet in a hard process. It carries the 
fraction x of the total energy E of the jet. Qq is the half opening angle of the jet while is the angle 
corresponding to the first splitting with energy fraction x <C z <C 1. 

2.1 Variables and kinematics 

The variables and kinematics of the process under consideration are the same as in section 3.1 of fill. 

2.2 Evolution equations for particle spectra at MLLA 

We define like in [ 1 ] the logarithmic partem densities 

Q(£) = xD Q (x), G(£) = xD G {x) 

for quark and gluon jets in terms of which the system of evolution equations for particle spectra at small 
x (see eqs. (42) and (43) of [ 1 ]) read 

Q(£, y) = S(£) + di' [ V dy'^il' + y') (l - -8{l' - £)) G(£', y'), (1) 

c J "0 

G(£, y) = 5(£) + J d£' J" dy'j^f + y') (l - aS(£' - £)) G(£' , y'), (2) 
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where 



4N r 



11 AT , 4 T A 



nt=3 



0.935. 



(3) 



The terms oc | in © and oc a in (O account for hard corrections to soft gluon multiplication, sub-leading 
g^qq splittings, strict angular ordering and energy conservation. 

2.3 Evolution equations; steepest descent evaluation 

The exact solution of © is demonstrated in [ 1 ] to be given by the Mellin's integral representation 



G(l,y) = (£- 



J J (2vri) 2 Jo u + s\{uo + s)u) 

{i+y+\)JJ 



V » „-As 



V + S 



V + s \u + s 



Ms) 



(4) 



where we have exponentiated the kernel (symmetrical in (u, v)) 



1 ln M^,_ As . 



/3(o; — i>) V^w + s) 



© will be estimated by the SD method. The value so of the saddle point, satisfying Sp- 
reads (see [7]) 



1 



v/^ + - v? - (w + f ) 



(5) 
= 0, 

(6) 



One makes a Taylor expansion of <x(s) nearby sq: 



<t(s) = (t(so) + ^"(s )(s - so) 2 + 0{(s- so) 2 ) , (t"(s ) = -!3\ 2 ^^ + {u-u) 



2 <0, (7) 



such that 



as /a;(^ + s)\ ' ^ Y i/ \ ' -As A>1 



Z/ + s\(cj + s) z/ 



^ + s 



= er(so) 



2 (f + so)vV'(so)| W*o 



a//3 



(8) 



The condition A » l=^a s /7r <C 1 3 guarantees, in particular, the convergence of the perturbative approach. 
Substituting © in (@) yields 



(2ttz) 2 (y + s )vV( s o)| V + so 



o/j9 



(9) 



where the argument of the exponential is 

4> (a;, z/, y) = u;^ + i/y + 



1 



ln^±^l-As . 



P{u) — v) (uj + s ) f 



(10) 



3 in 0, A appears to the power 3/2 > 1, which guarantees the fast convergence of the SD as A increases. 
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Once again, we perform the SD method to evaluate ©. The saddle point (ui , u ) satisfies the equations 
°"-> 1 ■lnr ( " + S : ) + / ,-a£±*£ = 0, (11a) 



dio /? (w — v) 2 ( w + s o) v (3lo(lu - v) (to - u) 



7r = y + 777 T2 ln 7 — i \ 15 — 7 \+ X l \ =0 - ( llb ) 



Adding and subtracting (II lab and dllbt gives respectively 



/5 + y + A) 



wo (^o + go) _ . wo + z/p + 2 s _ 
p" (w - v Q ) \ uj ' f J f3 (ujq - u f ~* (wo + so) ^o w - z^o 



y~ i = i57-. — — (— + — )-— ^ ln ~ i — -. ~ A — - — - — ; ( 12b ) 



(luq, vq) also satisfies (from ®) 

(wo + s ) (u + s ) = (13) 
One can substitute the expressions (11 lab and (II lbl) of £ and y into (flOl . which yields 

/ / /) \ 2 u; (fo + so) 

tp = <p{uj ,v ,l,y = — r-ln- ■ — . (14) 

(3 {U)Q - Vq) (lo + s ) vq 

Introducing the variables (n, v) [7| to parametrize (uq, vq) through 

coo (u ) = 1 (wo + so) (u + so) = -i^(^), (15) 

vT^+y+A) vpA 

one rewrites (fill) and d!2bt respectively in the form 

^, v ) = -LU e + y + X-V\) . . , (16) 

^/p V / smn |U — sinn v 



y — £ (sinh 2/z — 2/j.) — (sinh 2t> — 2u) 
y + ^ 2 (sinh 2 \i - sinh 2 t>) 

moreover, since ^0 — ^0 = (wo — so) — (^o — s o)> (/-*) v ) a l so satisfy 

sinh u sinh /z 



(17a) 



(17b) 



Performing a Taylor expansion of 4>{uj,v,£,y) around (ojq, z/q), which needs evaluating and 
(see appendix lA.lt . treating (Y + A) as a large parameter and making use of (fT5l provides the SD 

result 



G(£, y) « AA(/i, u, A) exp [— (V^ + y + A - a/A~) . r M " r + v - -^0* - u) 

Lyp V / sinn/j-sinni; p 

where 



(18) 



tf(p,v,\) = -(^+y+A) 



1 . I?) V4 / A W 2/3 



2 ^/TTcosh^De^/x,^) V^ + 2/ + A 
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with (see details in appendix IA.1I ) 



DetA{n,v) =(3(£+y+\f 



(fi—v) cosh [i cosh v + cosh fi sinh v — sinh \i cosh v 
sinh 3 /x cosh v 



(19) 



2.3.1 Shape of the spectrum given in eq. ffsl 

We normalize (fT8b by the MLLA mean multiplicity inside one jet [ 8 1 



n(Y) 



a»i 1 (Y + A 



A 



1 a 1 

~2/3 + 4 



exp 



2 

7^ 



The normalized expression for the single inclusive distribution as a function of 
ingly obtained by setting y = Y — I in (fT8b 



ln(l/x) is accord- 



ed, y) / /5 1 /2(y+A)3/2 



n(y) V ir cosh. vDetA(fi,v) 
One can explicitly verify that d20b preserves the position of the maximum l8lll9l H10I at 



Vy+\-\/x 



l)+ a ( ) 
sinh fx— sinh v J (3 



.(20) 



(21) 



as well as the gaussian shape of the distribution around d2TT > (see appendix IA.2t 



G(£,Y) 
n(Y) 



1/2 



tt^ [{Y + A) 3 / 2 - A 3 / 2 ] 



exp 



{£■ £-rnax ) 



(y + A) 3 / 2 - A 3 / 2 



(22) 



In Fig^we compare for y = 10 and A = 2.5 the MLLA curve with DLA (by setting a = in (J^UJi). 
The general features of the MLLA curve d20T > at A / are in good agreement with those of 1 2 1. 



0.2 



0.15 



MLLA, Y=1 0, Ji=2.5 
DLA, Y=10, \=2.5 




0.05 



9 10 



Figure 1: SD normalized spectrum: DLA (blue), MLLA (green); Y = 10.0, A = 2.5. 



The shape of the single inclusive spectrum given by d20b can easily be proved to be "infrared stable" (it 
has indeed a final limit when A — > 0). 
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2.4 Logarithmic derivatives 



Their calculation is important since they appear in the expressions of 2-particle correlations. 

Exponentiating the (I, y) dependence of the factor M in ( fflfl . we decompose the whole expression in 
two pieces 

V> = <p + Stp, (23) 
where ip, given in (U6t . is the DLA term for the shape of the distribution [7 1, and 

w = -U 1 + f) ln (^ + y + A ) - -f + f 1 + ^) u + 5 MQ(m, ^ ( 24 ) 

is the sub-leading contribution (in the sense that its derivative gives the MLLA correction), where 

Of v ) = P(£ + y + X ) 3 = sinh 3 (i 

' cosh v D et A(fjL,v) (fi — v) cosh [i cosh f+cosh/i sinh v — sinh// cosh v 

By the definition of the saddle point, the derivatives of (fT6t over I and y respectively read: 



(p e = U } = joe^, 
We introduce (see appendix IA.3t 



<Py = "0 = 70 



e 



a 1 d 

v ) = - jj + L (v, v ), L(fi, v) = -—]n[Q(n, v)], 



C(ii,v) = l + - + K{n,v) 



and make use of 



dv ( dfi 1 2 \ 

-=tanh^(coth^---/? 7o ), 



1 5 



-=tanh^ CO th / ,---/3 7o j, 



that follows from d!7bt . to write Sipg, 5tp y in terms of ^ 



= — 2 f ^ an ^ v ^(a*! u ) ) /^7o + ( £'/'• '') + ,au '' ''('<»th ///v'(//. r) 



9/i 



(25) 



(26) 



(27a) 



5ip v = — ( 1 + — + tanht;/C(/x, u) ] (3~/n + ( £(u, v) + tanhucoth/t /C(it, u) ) — — . (27b) 



Using dl7at and (I17bl) we obtain 
dfi 1 



;/?7o 2 



l + e"Q(/i,u) 



9y " 2^ 7 ° 



l + e-^Q(/i,u) 



(28) 



where 



Q{n,v) 



cosh /i sinh /u cosh v — (p — v) cosh t> — sinh v 



(fj, — v) cosh fj, cosh i; + cosh [i sinh u — sinh [i cosh t; ' 
which we have displayed in Fig|2](useful for correlations). 



(29) 
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Y+X=7.5 




l=ln(1/x) 

Figure 2: Behavior of Q(fi, v) as a function of £ = ln(l/x). 



Inserting (I27bt and d28l into (I27at gives the SD logarithmic derivatives of the single inclusive distribution 



ipi(fi,v) =7oe /i + -a7o e^Q(fj>,v)—tanhv — ta.nhvcothfi[l+e IJ 'Q(fi,v] 



- 2 H 



l+tanhu [1+K(p,v) )+C(/i,u) (l+e'*Q(/i,i;)) + 0(7o), (30a) 



ipy(/jj,v) = joe M — -O7o 2+e t>)+tanht> — tanhf coth^(l+e ^Q{^i,v 

-\pil\WBX&v(l+K(ix,v) \ -C(v,v)(l+e-»Q(fi,v))] + O( 7o 2 ) (30b) 



where we have introduced (L and K have been written in d67l and d68l ) 

C(fi, v) = L(fi, v) + tanhi; coth^l + K((j,, v)J . 
C does not diverge when /i ~ v — > 0. One has indeed 



(31) 



lim [L(fi,v) + tanhwcoth/iK( / u,u)] = lim - 



-H = 



as well as 



lim tanht; coth /U ( 1 + e ±tJ- Q([i, v 



lim 



q d 



l 



Y+X 



3/2 ' 



In (I30at and d30bl) it is easy to keep trace of leading and sub-leading contributions. The first 0(70) term 
is DLA [7] while the second (oc a — > "hard corrections") and third (oc /3 — » "running coupling effects") 
terms are MLLA corrections (0(7q)), of relative order 0(70) with respect to the leading one. In Fig|3] 
we plot (I30at (left) and d30bl> (right) for two different values of A; one observes that ip£ (ip y ) decreases 
(increases) when A increases. 

For further use in correlations, the logarithmic derivatives have the important property that they do not 
depend on the normalization but only on the shape of the single inclusive distribution. 
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1234567 0123456 



l=ln(1/x) l=ln(1/x) 

Figure 3: SD logarithmic derivatives ipg and ip y of the inclusive spectrum at Y = 7.5, for A = 1.5 and 
A = 3.5. 



2.4.1 "Limiting spectrum": A — > (Q = Aq C d) 

Since the logarithmic derivatives are "infrared stable" (see above), we can take the limit A — > in 
d3QaM30bl> 4 , and compare their shapes with the ones obtained in 0; this is done in Figs. |4] and |5l at 
LEP-I energy (EQ = 91.2 GeV, Y = 5.2) and at the unrealistic value Y = 15. 

The agreement between the SD and the exact logarithmic derivatives is seen to be quite good. The 
small deviations (< 20%) that can be observed at large i (the domain we deal with) arise from NMLLA 
corrections that one does not control in the exact solution. The agreement gets better and better as the 
energy increases. 




I=ln(1/x) l=ln(1/x) 

Figure 4: SD logarithmic derivatives tf}£ (left) and ip y (right) compared with the ones of Q] at Y = 5.2. 

It is checked in appendix 1A.41 that (fT8b satisfies the evolution equation 0; the SD logarithmic deriva- 
tives d30at and d30bt can therefore be used in the approximate calculation of 2-particle correlations at 
A ^ 0. This is what is done in the next section. 

3 2-PARTICLE CORRELATIONS INSIDE ONE JET AT A / (Q # A QCD ) 

We study the correlation between 2-particles inside one jet of half opening angle within the MLLA 
accuracy. They have fixed energies x\ = u\/E, X2 = Wz/E (uj\ > L02) and are emitted at arbitrary 

4 For this purpose, I17ai has been numerically inverted. 
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0.5 





>)/|, Exact Solution, Y=15 
Vl, SD, Y = 15, X°=0 



6 8 10 

l=ln(1/x) 



12 14 



3 p 

2.5 

2 - 

1.5 - 

1 - 

0.5 - 

_ 




V)/ y , Exact Solution, Y=15 
if SD,Y = 15,1=0 



6 8 10 

l=ln(1/x) 



12 14 



Figure 5: SD logarithmic derivatives ipg (left) and ip y (right) compared with the ones of JT| at Y = 15. 

angles 0i, 02- The constrain @i > 02 follows from the angular ordering in the cascading process. One 
has > ©i (see Fig. 1 of O). 

3.1 Variables and kinematics 

The variables and kinematics of the cascading process are defined like in section 3.2 of (Q- 



3.2 MLLA evolution equations for correlations 

The system of integral evolution equations for the quark and gluon jets two-particle correlation reads 
(see eqs. (65) and (66) of Q) 

Q {2) (£i,y 2 ,v)-Qi(h,yi)Q2(h,y2) = 1 f dt dy 1 2 (£ + y)\l--5(£-e l )\G^(£,y,r ] ), (32) 

c J J 



rh ry2 r -, 

G ( - 2) (e 1 ,y 2 ,v) ~ G 1 (£ 1 ,y 1 )G 2 (£ 2 ,y 2 ) = d£ dy 7 2 (£ + y)l-aS(£-£ 1 ) G^(£,y, V ) 

Jo Jo L J 



U2 



+ (a-b) dy^(£ 1 +y)G(£ 1 ,y + ri)G(£ 1 + V ,y). (33) 
Jo 



a is defined in (TJ} while 



AN r 



11 



^ c --n f T R [l-2^ 



2n 



rit =3 



0.915. 



(34) 



3.3 MLLA solution at A ^ 

The quark and gluon jet correlators C q and C g have been exactly determined for any A in Q by respec- 
tively setting Q( 2 ) = C q Q\Q 2 and G^ = C g G\G 2 into ( T3*2l and d3*3l . In the present work we limit 
ourselves to the exact MLLA solution which consists in neglecting all O(7o) corrections in equations 
(64) and (84) of fl]. 
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3.3.1 Gluonjet 

At MLLA, the logarithmic derivatives of tp d23l can be truncated to the saddle point derivatives tpi, tp y 
of (QUl- The MLLA solution of fTQ then reads (see (77) in Q) 



Cg-1 » x , A , r — i (35) 



MLLA 1 — 6 (y2l,£ + <P2,l) ~ $1 

1 + A + A' + Si 

where we introduce 



A = 7o 2 \V>l,e'P2,v + yi,y<P2,* J , (36) 

A' = 7o~ 2 (<Plj8lp2,y + tyl,y*P2,£ + Hl,W2,y + <-Pl,y5i>2,tj \ (37) 

X = In 1 + t , X^=-^r> Xj/ = -^-; (38) 
V l + Ay Xdy 



d36l is obtained by using d25l : 



(39) 



A(/ii,/i 2 ) = 2cosh(/ii-/i 2 ) = 0(1), (40) 
which is the DLA contribution Q, while d37l (see appendixlBji is obtained by using d25l . (I27at and (I27bt 

A'(/ui, At 2 ) = ; ; = O(7o); (41) 

7o 

it is a next-to-leading (MLLA) correction. To get d38l . we first use d40l . which gives 

l + 2cosh(//i-/i 2 ) V ^ ^ l + 2cosh(/xi-^ 2 ) V 0y ^ / ' 
and then d28t to get 

tanh^^ e^Qi-e^Q 2 tanh e ~^Qi - e~^Q 2 

X«-P7o 1 + 2cos h(^ 1 -/x 2 ) 2 ' Xy ~ ^ 7o l+2cosh(/x 1 -^ 2 ) 2 

which are 0(jq). They should be plugged into d39l together with d25b . which gives 

2sinh 2 ( ^~^ 2 ) f ~ ~ \ 

Jl = /?70 8 1 i^~(k^ \ Q ^ UVl)+Q ^ 2 ' V2) ) (43) 

it is also a MLLA term. For Q 3> Qq > Aqc d we finally get, 

1 + 2cosh(^i - fi 2 ) + A'(/ii,/x 2 ) + St 

where the expression for A' (1741 is written in appendix |E] It is important to notice that 5i ~ near 
£i f« ^ 2 (fii sa /i 2 ) while it is positive and increases as rj gets larger (see d29t and Fig|2]); it makes the 
correlation function narrower in \£\ — ^ 2 |. 
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3.3.2 Quark jet 

The MLLA solution of (|32j reads (see (93) in (TJ) 



C q -1 MLLA N c 



l + (6-a)(^ + 02,/) 



1 + A 

2 + A 



(45) 



Inserting d36l-(l39t into d45t we get 



c q (ii,h,y,\) 



MLLA 



N c ( 

1+ t c 



(Cg(e 1 ,£ 2 ,Y,X)-l 



1 + 



a)7o(e w + e 



/'2 



^ 1 + 2cosh( / ui - fi 2 ) 
2 + 2 cosh(^i — /i2) 



which finally reduces (for Q » Qq > Aqcd) to 



MLLA 



n — - 



\ 1 _|_ P M2 

C g (4,4T,A)-l)+-(6-a)7o r ' 



+ cosh(/ii - /x 2 ) 



(46) 



3.4 Sensitivity of the quark and gluon jets correlators to the value of A 

Increasing A translates into taking the limits /3, Aqcd —>0(Y = £ + y<^\,Q^>Qo^> Aqcd) in 
the definition of the anomalous dimension via the running coupling constant (70 = 70 (a s ), see (44) in 
[ 1 1). It allows to neglect £, y with respect to A as follows 



1 



1 



£+y<.\ 2 

W + y + A) ~ lo = J\ 



(47) 



such that 70 can be taken as a constant. Estimating © in the region A ^> 1 44> s <C 1 needs evaluating 
the kernel 



1 f uj (u + s) 



v + s \ (u + s) V 



V 



V + s 



a/0 „ 1 / \ l//3(w-i>) 

<<11/ l+^ s ) [1 



1 + 



1 1 



V \ UJ 



s 11/1 
"^ + 2!^U- a 



2 s 2 



UJV 

11/1 



V / 



f3 2 + 3! 1/ 3 \uj a ) /3 3 + 



3 S 3 



Integrating ( 148 1 over s, using d47t and J c 



00 „n „-As _ n 



S e 



A 



tt, we get 



\ + -(--a]^T + ^(--a] (-M +4(--o| f^r 



V \ UJ 

1 



/3A ^ 2 



/3A 



z^ 3 \ u 



/3A 



^-7o (1/w-o)' 
which, after inverting the Mellin's representation (132) of [ 1 1, gives 

G(l, y) s ~ exp(27o y^y - a7o2/)- 
Taking the same limit in dl7at and d!7bt gives respectively 

z+y<.\ 1 y — I 



y-i i+y<\ 1 y 
tanh u =^ u = - m — , 



jl — V 



2 A 



~ -u. 



Furthermore, we use d50l to show how d23b reduces to the exponent in d49l 5 

2 — 1> l+y<^\ 



y/~j3 sinh /i — sinh v 



270 v^y, 



5 we set (3 = in <30a> . I30bl and only consider terms cx a 



(48) 



(49) 



(50) 
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VQ 



a/p 



= -~ln[l 

2(3 

t+y<.\ 2 
ps -a7 y. 



Thus, since /j, = \ In ¥ d50b . d30al) and (I30bt simplify to 



1 a £ + y a , 
2/3~ /T^ 



(51) 



7o e^ = 7 



% « Toe ^ - a7 = 7 01 



(52) 



Therefore, taking the limit (3, Aqcd —> (A — > oo) leads to the simplified model described in section 
4.2 of IT1. Setting, for the sake of simplicity, l\ ps £2 in (l44ld45l . where 5\ vanishes, we obtain, in the 
high energy limit 



C g (£,y)~l + -, 



where 



1-2 [b--a]i/>t(£,y) 



, c q (i,y)-i+ 



Nc 

c F 



1 1 /5 



3 2 13 



a+b 4>i(£,y) 



1 1 / T R T R C F T R Cl\n f =3 



0.6, 



(53) 



5 L 2 
-a + 6 = - 

3 9 



Tr TrCp T R Cl\ n/=3 



(54) 



Thus, when A increases by decreasing Aqcd, ipi oc 70 decreases and the correlators d53t increase. For 
LHC, a typical value is Y = 7.5 and we compare in Fig. |6] at fixed Qq, the limiting case A pa 
(Qo ~ Aqcd ~ 253MeV) with A ps 1.0 (Aqcd = 100 MeV) and A « 2.3 (A QCD = 25MeV). As 
predicted by d53l . the correlation increases when Aqcd —> at fixed Qq. 

It is also sensitive to the value of Qq. As seen in d53l . since y = In M — £, if one increases Qq (since 
Aqcd is fixed, 70 does not change), thereby reducing the available phase space, the correlators increase. 
This dependence of the correlators at fixed Aqcd is displayed in FigQfor 0.3 GeV < Qq < 1.0 GeV at 
£ 1= £ 2 = 3.0 (soft parton). 





Figure 6: Varying A at fixed Qq; Aqcd dependence of C g (left) and C q (ri| 



;ht) 



In the simplified model which leads to d53l . C g and C g respectively go to the asymptotic values 4/3 and 
1 + N c /3Cf- This is however not the case in the general situation (3 7^ 0, as can be easily checked by 
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1.35 
1.3 

1.25 
1.2 

a 115 
if 1-1 

1.05 
1 

0.95 
0.9 




Figure 7: Varying A at fixed Aqcd = 253 MeV; Qo-dependence of C g (ieft) and C q (right) at £\ = £2 
3.0 



using (I30at and (I30bt : for example, near the maximum of the distribution (/i ~ v — > 0), a contribution 
oc A 3 / 2 /[(y + A) 3 / 2 — A 3//2 ] occurs in the term proportional to in d53t that yields negative values of ipi 
when A increases. 



3.5 Extension of the Fong and Webber expansion; its limit A = 

In the Fong-Webber regime, the energies of the two registered particles stay very close to the peak of the 
inclusive hump-backed distribution that is, \£i - £ max \ < a oc [(Y + A) 3/2 - A 3/2 ] 1/2 (see d2~2t). 

Near the maximum of the single inclusive distribution l\ ~ £2 ^ Y/2 (/u, v — > 0, see appendix IA.2t 



A \ 1/2 
lim C = I — — 

^-►0 V Y + A , 



v r, 3 «i i- 2 1 / A x 3/2 
lim Ki = 5—? — 5-, km Q = - + - , 



where C, ifj and Q are defined in (Bit . (I68t and (I29t . Keeping only the terms linear in /j and the term 
quadratic in the difference (/xi — ^2), one has 



A + a' ~ 2 + {m- 112? - 070 (2 + m + m) - /?7o 



2 + 3 



A 3/2 



(y + A)3/2 _ A3/2 



(55) 



and 



«i~* 2 ~y/2 1 2 
<h ^ t:P7o(mi-M2) 



2 + 



A 



y + A 



3/2' 



(56) 



<5i can be neglected, since 7o(mi ~~ M2) 2 *C (/ii — /i2) 2 *C 1. Then, in the same limit, (PHI . ( Pol l become 



1 - 670(2 + m + /i 2 ) 



3 + — /J2) 2 - a7o (2 + Ml + M2) - /?7o 



2 + 3 



A 3 / 2 



(y + A) 3 / 2 - A 3 / 2 



(57) 



C° q (£ 1 ,£ 2 ,Y,\) * 1 + ^ 



C°(^,4,y,A)-l)+i(6 - o), n (2 + //L+//0) 



(58) 
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Using ( E)4T ) one has 

(/^i-/"2) 2 ^ 9 



Y + X 



[(Y + A) 3 / 2 - A 3 / 2 ] 
such that the expansion of ( 1571 . d5%t in 70 oc ,Va7 reads 



(y + a) 3 / 2 - a 3 / 2 



C g Q (^ 2! y,A) 



V + A) 1/2 (li-l 2 ) N 

(y + A) 3 / 2 - a 3 / 2 



h f 2 + (y + A) 1/2 y V 1 



+ 



1 / 2 



+ 



A 3 / 2 



3 I 3 (y + A) 3 / 2 - A 3 / 2 



/37o + ( b - -a 



3 (y + A) 3 / 2 - A 3 / 2 y V3 

(y + A) 1/2 (i x + i 2 ) x 
(y + A) 3 / 2 - a 3 / 2 



a - 61 70 
7o + 0(7o), (59) 



c ? Vi,4,y,A)^i+^+^ 

A 3 / 2 



_ / (y + A) 1 / 



(y + A) 3 / 2 - a 3 / 2 




+ 3 I 3 + (y + A) 3 / 2 - A 3 / 2 J ^° + 4V3° + 6 ] 



(y + A) 3 / 2 

(y + A) 1 / 2 (ii + i 2 A 
(y + A) 3 / 2 - a 3 / 2 




-a+M 70 



(60) 



Therefore, near the hump of the single inclusive distribution, (l4*4l .(l4*6Tl behave as a linear functions of the 
sum (£1 + £2) and as a quadratic functions of the difference {£\ — £2). At the limit A = 0, one recovers 
the Fong-Webber expression (3). 



3.6 Comparison with the exact solution of the evolution equations: A = 

In Figsj8]we compare the SD evaluation of the gluon correlator with the exact solution of [ 1 1 at A = 0. 
The difference comes from sub-leading corrections of order 7q that are not present in (PPfli . For example, 
— /?7 2 ps —0.2 at y = 5.2 occurring in the exact solution (69) of {Q is not negligible but is absent in 
d44t and That is why, the SD MLLA curve lies slightly above the one of Q at small £\ + £ 2 . 

The mismatch becomes smaller at Y = 7.5, since — /?7 2 ~ —0.13. However, when £\ + £2 increases, 
the solution of [ 1 1 takes over, which can be explained by comparing the behavior of the SD MLLA $i 
obtained in d43t and 6 C , 5 C in Namely, while 5\ remains positive and negligible for £\ £2, 5 C , 5 C 
decrease and get negative when £\ + £2 — > 1Y, see Fig|9](left), which makes the correlations slightly 
bigger in this region. As \£\ — £ 2 | increases, 5\ is seen in Figj9] (right) to play the same role as 5 C , 5 C 
do in the solution [ 1 1 and therefore, to decrease the correlation. The agreement between both methods 
improves as the energy scale increases. A similar behavior holds for the quark correlator. 

In Q, strong cancellations between the MLLA S± and the NMLLA 62 were seen to take place, giving 
very small 5 C and 5 C ; this eased the convergence of the iterative method but raised questions concern- 
ing the relative size of MLLA and NMLLA corrections and the validity of the perturbative expansion. 
However, since 5i is itself, there, entangled with some NMLLA corrections, no definitive conclusions 
could be drawn. The present work and Fig. |9j by showing that, below, 5 C and <L of lITT play the same 
role as the pure MLLA 5\ which is now calculated, suggests (though it is not a demonstration) that the 
perturbative series is safe. It is indeed compatible with the following scheme: in Q, the pure MLLA part 
of 8\ is the same as that in the present work; the cancellations in Q occur between NMLLA corrections 
included in 5\ and 62 ; these are eventually of the same order of magnitude as MLLA terms, but they are 
only parts of all NMLLA corrections; this leaves the possibility that the sum of all NMLLA corrections 
to 5\ and all NMLLA terms of <5 2 are separately smaller than the pure MLLA terms of S\, that is that 
strong cancellations occur between NMLLA corrections, the ones included, because of the logic of the 
calculation, in [ 1 1, and those which were not be taken into account. 
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Y=5.2, X=0 



Y=7.5, X=0 




Figure 8: Comparison between correlators given by SDandin[l],atA = 0. 



3.7 Comparison with Fong- Webber and LEP-I data; how A = is favored 

Let us consider, at the Z° peak Y = 5.2 (EQ = 91.2 GeV at LEP-I energy), the process e + e~ — > qq. As 
can be induced from FigJS] the results obtained in the present work by the (approximate) SD method are 
very close to the ones obtained in subsection 6.5 of [1 1 by the exact solution of the evolution equations. 
Accordingly, the same comparison as in [ 1 1 holds with respect to both Fong & Webber's results [ 3 1 and 
OPAL data 0. 

It is also noticeable that, since, at A = 0, correlations already lye above (present) experimental curves, 
and since an increase of A tends to increase the predictions, the limiting spectrum stays the best candidate 
to bring agreement with experiments. 



4 CONCLUSION 



Let us, in a few words, summarize the achievements, but also the limitations of the two methods that 
have been used respectively in Q (exact solution of MLLA evolution equations) and in the present work 
(steepest descent approximate evaluation of their solutions). 

Achievements are threefold: 

- in Q, MLLA evolution equations for 2-particle correlations have been deduced at small x and at any 
A; their (iterative) solution can unfortunately only be expressed analytically at the limit A — > 0; 

- by the steepest descent method, which is an approximate method, analytical expressions for the spec- 
trum could instead be obtained for A ^ 0, which enabled to calculate the correlation at the same level of 
generality; 

- one could move away from the peak of the inclusive distribution. 
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So doing, the limitations of the work of Fong & Webber have vanished. Their results have been recovered 
at the appropriate limits. 

The two methods numerically agree remarkably well, despite an unavoidable entanglement of MLLA + 
some NMLLA corrections in the first one. 

The limitations are the following: 

- the uncontrollable increase of a s when one goes to smaller and smaller transverse momenta: improve- 
ments in this directions mainly concern the inclusion of non-perturbative contributions; 

- departure from the limiting spectrum: it cannot of course appear as a limitation, but we have seen that 
increasing the value of A, by increasing the correlations, does not bring better agreement with present 
data; it confirms thus, at present, that the limiting spectrum is the best possibility; 

- the LPHD hypothesis: it works surprisingly well for inclusive distributions; only forthcoming data will 
assert whether its validity decreases when one studies less inclusive processes (like correlations); 

- last, the limitation to small x: it is still quite drastic; departing from this limit most probably lye in the 
art of numerical calculations, which makes part of forthcoming projects. 

Expectations rest on experimental data, which are being collected at the Tevatron, and which will be at 
LHC. The higher the energy, the safer perturbative QCD is, and the better the agreement should be with 
our predictions. The remaining disagreement (but much smaller than Fong-Webber's) between predic- 
tions and LEP-1 results for 2-particle correlations stands as an open question concerning the validity of 
the LPHD hypothesis for these observables which are not "so" inclusive as the distributions studied in 
0. The eventual necessity to include NMLLA corrections can only be decided when new data appear. 
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A DOUBLE DERIVATIVES AND DETERMINANT 

A.l Demonstration of eq. dl9t 

We conveniently rewrite (11 lab and il lbt in the form 

86 2uj — v 



8uj u — v uj — v uj — v 
86 uj — 2v 



1 



-y 



+ 



dv uj — v uj — v uj — v 
The Taylor expansion of dlOt in ^ reads 



+ A 



u — v (3uj(uj — v) ' 

uj + 2sq 1 
uj — v [3u(uj — v) 



(61) 
(62) 



6(uj,v,£,y) 6(uj ,u ,£,y) + ^ ^-^ (w , Vo)(w - uj ) 2 + - (uj ,vq)(u - v Q f 



1&< 
2 du 2 



2 



+ duidv 1/0 ^ ^ ~ ~ ^ 

The expressions of the second derivatives follow directly from (161 1 and d62l 



eta 2 (w — i/) 2 



(^+y+A) + 



2w — ^ 



+ 



(w-^) 2 (3uj 2 (uj - vf f3(uj-u) 2 (2s + uj + u) 



(63) 



2 



du 2 (uj — v 



uj , a; — 2z/ 

1 +y+A) + t 777 + ~ — ~ + 



(uj-u) 2 (3u 2 (uj-u) 2 P(uj - v) 2 (2s + uj + v)' 



2 



ii+y+X) 



+ 



1 



dujdv (uj-v) 2 ' y ' (uj-v) 2 [5uj(uj-v) 2 (3(uj-v) 2 (2s + lu + v)' 
^ and its solution can be written in the form 

2vr 



G~ / / d z ve-2' 



VDetA 



where 



uj 



v = (uj, v), v 1 



DetA = Det 



/ d 2 (f> a 2 <t> \ 

du; 2 dujdu 

d 2 <t> d 2 4> 

\ dudu! du' 2 I 



cP^d 2 ^ _ f d 2 6 x 2 
duj 2 dv 2 \dujdv 



An explicit calculation gives 

DetA = (£ + y + \f 
which, by using (Tf5l leads to ( Tf9t . 



0(u + v)6-A 



+ 



4(cj + v) 



(uj-v) 2 (uj - v) 2 (2s + uj + v) 
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A.2 DetA (see eq. dl9b ) around the maximum 

This is an addendum to subsection 12.31 

imax written in (ED is close to the DLA value Y/2 [ 7 ] [ 8 ] [ 9 ] . We then have /j~t;^0for^y~ F/2. 
In this limit, d!7al) and dl7bl) respectively translate into 

„ m,w-o 2 (y + A) 3/2 - A 3 / 2 



while 



9£ " (y + A)3/2 _ as/2 
should be used to get (122b . An explicit calculation gives 



% = _3 7 iY + f\„ (65, 



lim 



//?i/2(y + A)3/2 



1/2 



/v>->oy TtDetA(fi,v) y 7T ^[(Y + A) 3 / 2 -A 3 / 2 ] 
where 



/?(y+A) 3 
1 



fci-t,) (l + i/, 2 ) (l + ^ 2 ) + (l + i/, 2 ) (t^ 3 )-^ 3 ) (l + ^ 2 ) 



-/?(y+A) 3 i 



1 



/?(y + A) 2 



A \ 3 / 2 ' 



(66) 



A.3 The functions L(ii, v), K(n, v) in eq. d26b 

An explicit calculation gives 

3 cosh 1 (p, — v) cosh f sinh /i + sinh v sinh 



and 



L{fi,v) 



2 sinh fi 2 (fi — v) cosh /j, cosh v + cosh /j, sinh v — sinh cosh v ' 



sinh f 



(/i — v ) cosh /i — sinh ^ 



(/x — f) cosh /u cosh t> + cosh [i sinh t> — sinh [i cosh t> 



(67) 



(68) 



A.4 A consistency check 

Let us verify that the evolution equation © is satisfied by d20t within the MLLA accuracy. Differentiat- 
ing © with respect to I, y yields the equivalent differential equation 



G iy = 7 2 (G-aG £ ) + O(^G) 



that can be rewritten in the form 



i>li>y + ^ly = To C 1 - a ^t) + <9(7o) ; 



(69) 



we have neglected next-to-MLLA corrections O(7o) (of relative order 7 2 ) coming from differentiating 
the coupling 7 2 in the sub-leading ("hard correction") term oc a. 
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We have to make sure that d^9l holds including the terms 0(jq). In the sub-leading terms we can set 
ifr — > (p (see (l2"5t): 



(ipt + 5ipj>)(ipy + 5ip y ) + ip £y = 7^(1 - a</^). 
Isolating correction terms and casting them all on the l.h.s. of the equation we get 



(70) 



aioW + Wltyy + Vytyl ]+<P£y = To _ WW 



(71) 



By the definition d25l of the saddle point we conclude that the r.h.s. of dTTT i is zero such that we have 



that is, 



dujQ 



dy 



0. 



First, we select the terms oc a: 



(72) 



(73) 



"To 



1-1 1 1 

- — Q — — tanh v + — tanh v coth \i + — tanh t> coth Q 



+—Q tanhfe tanh v coth u e ^ tanh t> coth uQ 

2 2 2 2 

= [— tanh v cosh + tanh v coth sinh /j] = 0. 
From (fT3l one deduces 

= 2^ 7 ° g ' 

that is inserted in d73l such that, for terms oc /3, we have 



-/3To 3 



-e^ + - tanhu(l+iv")e^ - -C ef* - -CQ + -e^ + - tanhy( 1 + A' ) < ->' 



-/3To 3 



(•( >sl) // + l.<n i h rc-( )s] i // j 1 + K ) — C sinh // — -Q 



which gives 



"/?To 3 

Constructing (see (129b and appendix lA.3t 



cosh /i — sinh (j,L — —Q 



Q(fi, v) — 2 cosh = —3 cosh /i + sinh /i 



(fj, — v) cosh t) sinh fj, + sinh i; sinh 
{fj, — v) cosh /i cosh v + cosh /i sinh v — sinh /i cosh v 



-2 sinh iiL(/j,, v) 



we have 



-/3To 3 



cosh /x — sinh /i L — — Q 



= 0. 
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B ANALYTICAL EXPRESSION OF A (^i, /x 2 ) OBTAINED FROM EQ. (EJ) 

Replacing d30at(l30bt in PTt and neglecting terms of relative order 0(7q) which are beyond the MLLA 
accuracy, we obtain 



A 



-a-yo 



7o 

e Mi _|_ e M2 _ s i n h(^ 1 — /i 2 )(Qi — Q2) + cosh hi tanh 1*2 + cosh ^2 tanhui 
— sinh [i\ tanh V2 coth H2 — sinh ^2 tanh v\ coth Hi 
+ sinh(^i — ^2) ^ tanh f 1 coth fi\Qi — tanh t>2 coth 

cosh/ii - sinh^iC 2 ) + ( cosh/i 2 - sinh/^Ci ) + sinh(/ii - ^XCiQi - C2Q2 



+ cosh fj>i tanh ( 1 + K2 ) + cosh /i2 tanh v\ ( 1 + K 



(74) 
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